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THEORETICAL MOTIONS OF HYDROFOIL SYSTEMS 
By Frederick H. Imlay 

SUMMARY 


Resul'bs are presented of an investigation that has "been under- 
taken to develop theoretical methods of treating the motions of 
hydrofoil systems and to determine some of the Important parameters. 
Variations of parameters include three distrlhutions of area 
"between the hydrofoils, two rates of change of downwash angle with 
angle of attack, three depths of immersion, two dihedral angles, 
two rates of change of lift with immersion, three longitudinal 
hydrofoil spacings, two radii of gyration in pitching, and various 
horizontal and vertical locations of the center of gravity. Graphs 
are presented to show locations of the center of gravity for stable 
motion, values of the stability roots, and motions following the 
sudden application of a vertical force or a pitching moment to the 
hydrofoil system for numerous sets of values of the parameters. 

The lateral stability of tandem-hydrofoil systems is briefly 
discussed, and values of the latersil stability roots are presented 
foi- two longitudinal hydrofoil spacings and two vertical locations 
of the center of gravity. 

The analysis indicates that if only the longitudinal motions 
of a hydrofoil system are of interest the present theory should 
provide satisfactory predictions. M adeqmte theory for the 
lateral motions, however, must treat the longitudinal and lateral 
motions in combination. The conclusions based on the investigation 
are that a large longitudinal spacing bet^^en the hydrofoils, a 
large rate of change of lift with depth of immersion, and a horizontal 
location of the center of gravity near the center of the region of 
stable locations are important contributions in the attainment of 
desirable characteristics for the longitudinal motion. An appendix 
gives an outline of the methods of theoretical treatment used and 
presents methods used in computing the required stability derivatives. 
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INTRODUCTION 


The use of hydrofoils as an alternative to planing "bottoms or 
hulls for the support of craft operating on the surface of water has 
been of interest for some time. (See reference 1 .) Guidoni 
advocated the use of hydrofoils as a means of improving the talce- 
off and. rough-water performance of seaplanes ns early as I9II. 

(See refei'ence 2 .) Some of the advantages claimed for hydrofoils 
over planing "bottoms are a "better ratio of lift to drag on the 
water and less sensitivity to irregularities of the water surface . 

In addition, if hydrofoils are used, the hull lines can "be designed 
to favor good aerodynamic rather than good hydrodjnamlc charac- 
teristics, and "by retracting the hydrofoils the aerodjaiamic 
perfomance can "be even further improved. In spite of the evident 
advantages of these devices and the attention that they have 
received, no pu"blished work is laaovTi to exist on the sta"bility of 
motion for systems employing hydrofoils . 

The present paper deals theoretically with the "behavior of a 
system supported solely "by hydrofoils and is a first approach to 
the pro"blem of developing methods of theoretical treatment for the 
more general case where the interaction of hydrofoils, hull, and 
aerodjoiamic surfaces have to "be taken into account. The treatment 
is "based on the theory of small oscillations and involves assumptions 
customarily made in applying the theory. (See reference 3.) 

Definitions of all sjoabols used are listed at the beginning 
of the appendix. 


LONGITUDINAL MOTIONS 


The longitudinal motions of a number of hypothetical hydrofoil 
systems were investigated by raesns of calculations based on the 
theoretical treatment presented in the appendix. All the computa- 
tions were for systems composed of two similar h;,'drofoils of rectangular 
plan form and rectan,gular tips. The hydrofoils were arranged in 
t;mdem and had an aspect ratio of 6 and a total hydrofoil area of 
0.188 square foot. (See fig. 1 .) The systems were assumed to have 
a mass of 0 .256 slug and to operate at a velocity of 20 feet per 
second in water having a density of 1 .97 slugs per cubic foot . The 
mass of the system was assumed to include all items such as structure 
and additional mass effect. For systems with dihedral the hydrofoil 
area, aspect ratio, and span were based on the nart of the hydrofoil 
immersed during the initial vindlsturbed motion, although \inwetted 
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parts of the hydrofoils were asstmied to pro.lect above the water far 
enou,<th to ensrre that the tips were never immersed during distiirhed 
motions. (See fig. 2.) Most of these dimensional characteristics 
of the hj''drofoil systems vrere chosen to facilitate comparison of the 
theoretica]. motions with the results of contemplated experimental teats. 
Changes in the other parameters were made to determine their effects 
on the stable regions, the stability roots, and the motions resulting 
from disturbances . 


Effect of Parameters on Stable Pegions 

The stable region, as used in the pi*esent paper, indicates 
pe;nnissible locations of the center of gra.vity rela.tive to the 
hydrofoils if the longituidinal motions are to be stable. The stable 
region alone, however, gives no quantitative indication of the 
degree of stability. The stable region is bounded by lines that 
are the loci of center-of -gravity locations for which neutrally 
stable longitudinal motions occur. The positions of the boundary 
lines, and hence the size of the region, vary with changes in the 
parameters of the hydrofoil system and thus suggest voTiations of 
the parameters that may be of practical interest for more detailed 
study . 

The type of mstable motion occuiTing Just outside the boundaries 
has been noted for each of the stable regions in figures 3 to 9 ; 
thus, for each stable region, center-of-gravity locations beyond the 
rear boundary lead to an unstable divergence, and in most cases 
unstable oscillations occur for locations beyond the front boundary. 
The rear boundary is always located farther- to the rear of the front 
hydrofoil than would be the case for a similar pair of airfoils 
because of the additional damping Introduced as a result of the 
sensitivity of the hydrofoils to depth of immersion. 


In addition to the selection of a center-of-gravity location 
that lies within the stable region in order to meet the requirements 
for stability, certain supplementary practical factors must be 
considered. For example, negative lift on either hydrofoil shoiad 
be avoided; otherwise momentary imcovering of the hydrofoil (as by 
a wave trough) will be followed by nosing-over if the rear hydrofoil : .. 
is operating at negative lift, or nosing-up if the front hydrofoil 
is operating at negative lift. Fvirthermore, the longitudinal 
location of the center of gravity is also restricted by the maximum 
positive lift obtainable, and may be Infliienced by the desirability 
of operating the hydrofoils near their maximum lift-to-drag ratios. 

The net effect of such restrictions Is to reduce the usable part of 
some of the computed stable regions shown in figures 3 to 9 . 
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In the pi-’esent study, where the effects of power are neglected, 
the vertical center-of -gravity location selected appears to he of 
secondary importance, lew locations heing somewhat advantageous. 

The effects of -oover, however, will undoubtedly have an important 
hearing on the choice of the vertical center-of -gravity location. 

Dlstrihution of ar ea.- The effect of the distribution of area 
between the two h 3 'drofoils on the extent of the stable region is 
shown in figure 3. The plan-form arrangements assumed for the 
three distributions treated are shown in figrure 1. In arrangement 1 
the hydrofoils were identical; in the other two arrangements the 
ratio' of the distribution of area was l:k and the arrangements 
differed only in the location of the larger hydrofoil. All the 
arrangements had the same total hydrofoil area of 0.108 square foot. 
The horizontal distance between the assumed hydrodynamic centers of 
the hj'drofoils for all arrangements was 10 .OC]^, where ci is the 
chord for the arrangement with two equal hydrofoils, and the assumed 
hydrodynamic center was located at the quarter-chord point of the 
center section. All the hydrofoils were assumed to be immersed 1 .0C]_ 
at the hydrodynamic center during the initial undisturbed motion. 

figure 3 shows that the configuration with the small surface 
ahead ( ai’rengement 3) gave the largest useful stable region. The 
rearwa.rd extent of the stable region for the arrangement with two 
hydrofoils of equal area (ari’angement 1) was considered adequate, 
however, end becaiise this arrangement permitted certain simplifica- 
tions in the calcAilations, it was used for the rest of the work. 

The configuration having the main surface ahead (arrangement 2) 
would, from theoretical considerations, be the most efficient 
arrangement for developing lift but has a considerably more limited 
range of stable- center-of -gravity location than do the other 
arrangements. ■ . i ■ . 


Rate of change of do^mwash .- In a tandem- hydrofoil system the 
downward velocity produced In the fluid by the front hydrofoil reduce® 
the effective angle," of attack of the rear hydrofoil by the amount of 
the doniifash angle e. The dewmrash an.gle is a function of the lift 
On the front hydrofoil and hence varies with angle of attack. The 
rate of change of downwash angle with angle of attack, which is the 
factor cf Interest from the standnolnt of stability, will be represented 
by the smbol The value of will probably be intermediate 

'2 

between zero and the theoretical ultimate max.imum e = * 

’^'^-1 Sai 

but to determine the value accurately would require .an investigation 
of downwash near a free surface . Corresponding limiting values of e, 
which are given instead of in the figures for the sake of 

brevity, are zero and twice the Induced angle of attack ai . In order 
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to show the • inflnence . of the I'ate of chen^je . of ..downwash on the 
natiire of the stable re^^lon, compf.tatior.s were made for these two 
extremes, and the resuats for a system hawing? two equal hydrofoils 
are shovn In fl,gare 1 , An increase in the variation of domwash 
with a shifts both boundaries forxrard without appreciably altering 
the size of the stable region. 

The effect of dovmwach for the other hydrofoil arrangements was 
found to be similar to that indicated by flgi.i'.re ^ for the arrangement 
with two equal hydrofoils. Because there was no pronoxmced change 
in the size of the stable region with ch.ange in doTOwash, the 
condition of zero rate of change of downwash with a was assumed in 
most of the remaining calculations. 

The tiuTe boundaries of the stable region for the system treated 
in fi.gu.rs U lie somewhere within the bends, defined bv the boundaries 
.for e = 0 and e = 2a^, but' accurate definition of the boundaries 
requires that e be ImCnm. Conservative estimates will b© obtained, 
when the value of e is not Imown, if the assumptions are made that 
0 — for computing the locaticzi of the rear boundary and that 

6 = 0 for the front boundar;y . ■ ' ' , 

The lift .and drag obtained from a hydrofoil 
depend upon the depth of immersion Zq of the hydrofoil in the water. 
Because appreciable change in the depth of Immersion may occur ■under 
normal operating conditions, computations of the stable region were 
made for Immersion depths of 0 . 5 ci, l.Oci, and 1 . 5 g,i. (See fig. 5.) 
Limits of the stable region were not aJ.teied to anj^ Important extent 
by the assumed changes in the depth at which the hydrofoils operate. 

The effect on the stable region of Increasing 
the dihedral aagle P of the hydrofoils from 0 ® to 30° is oho'v'/n in 
■figure 6 , Both the front and the rear boundaries of the stable 
region were affected by the dihedral in such a way that the increase 
In dihedral increased the size of the stable region. 

Increasing the dihedral angle from'O® to 30° resulted in an 
associated increase in vertical damping. It appeared reasonable 
that the Improved stability obtained by changing the dlhedrcJ 'mlglit 
have resulted from this increased vertical damuln.gj consequently 
the effect of arbitrarily increasing the vertical ‘damning for the 
hy..rofoils with a dihedral angle of 0° was studied and the results’ 
are discussed In the next section. ' 

lift, with im me-sion.- If the depth of immersion 
0- a hyduofoii is changing, the lift is also chan-glng, and the rate 
of change- can be expressed by the vertical-damping derivative 5 Ct / 8z ' . 
It is believed that the increased stability which accompanied the 
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increase in dihedral angle from 0° to 30 ° (discussed in the preceding 
section) may have heen "brought about "by the resulting increase in 
the v8-lue of ' • Inasmuch as a fu?cther increase in dihedral 

angle would decrease the value of the derivative, an explanation of 
the increase in ' when the dihedral was changed from 0° to 3CP 

may "be of interest. 


In ord.er to avoid the mathematical difficixlties of treating 
discontinuous derivatives the assumption was made in the present 
studj’', for the case of hydrofoils with dihedral, that a normally 
inactive part of the hydrofoil extended sufficiently far above the 
water surface to keep the hydrofoil, from being completely immersed 
at eny time during distu.rbed motion. (See fig. 2.) As a result of 
this ass\maption, hydrofoi.ls with dihediaal have a larger variation of 
lift with change in depth of irar.iersion. than do hydrofoils with 
0° dihedral because of the increased, area brought into action when 
the hydrofoil sinks deeper into the water. This variation in active 
area becomes greater as the dihedral angle becomes smaller. 

The effect on the sta.ble region of increasing the value of 
oCl/o2 ' for each hydrofoil to twice the value that the hydrofoils 
had with 30 ° dihedral, but having other characteristics the same as 
for 0° dihedral, may be seen by comparing figures 6 and ?• Doubling 
the value of ^Cj^/bz ' shifts the rear boindary of the stabl.e region 
back considerably and produces nronounced changes in the fron"; 
boundary. The fonaer boundaj'v for lunstable oscillations now becomes 
an unstable "hump" in the region with a new front boundary ahead of 
the hump. The new forward boundary represents conditions for an 
unstable divergence, but the boundarj'’ is too far ahea,d of the front 
hyd?.'’ofoil to be of an 5 '’ practical interest. 

Longitudinal h.vd.rofoil spacing .- The effect on the stable 
region of increasing the longitud'inal spacing of the hydrofoils from 
lOci 20c]_ is shown in figure 8. The larger spacing results in a 
very large increase in the stable region and in the replacement of 
the front boundary that indicated unstable oscillations by a new 
front divergent boundar.y. The new front boimdary is well ahead of 
the front hydrofoil, which is the practical limit of forward center- 
of - gravity loc .a,tion . 

The absence of a boundary for oscillatory instability for the 
system with a spa,clng of 20c^ suggests that the large amoun.t of 
damping in 'oitching for this spacing, rel.ative to the pitching radius 
of gyration Ky, might result in overdamping and thus prevent the 
system from having any oscillatory motion. Calculations with Ky 
reduced to give a similai" relation between inertia and daraping for 
the small spacing of lOc-j^, made to check the -hypothesis, showed that 
oscillations were still b"bta.ined) thus, it appears that the absence 
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of 'unsta'ole oscillations for the lai’ger spacing does not signify 
ina'oility of the system to have transient oscillations. 

The pronounced increase in the size of the stable region when 
the longitudinal spacing of the hydrofoils is increased Indicates 
that a large spacing is desirable in order to minimize the effects of 
unavoidable changes in center-of- gravity location encountered in 
practice. In a previous section entitled "Distribution of Area" 
a spacing of 10c]_ was used in the calculations made to study the 
effects of distribution of area. If a larger snacing had been 
used, it would possibly have resulted in a sufficient gain in the 
size of the stable region for the arrangement with the large 
hydrofoil forward to make this configuration of practical value. 

gadius o f gyration in pitching .- The marked increase, in 
permissible horizontal center-of- gi’avity movement when Ky is 
redxiced is indicated in figirre 9> whei”e the stable range of 
horizontal center-of -gravity location is shown for zero vertical 
elevation of the center of gravity with Ky reduced to one -fourth 
the value used previously. The pronounced effects of reducing Ky 
indicate that increased values of Ky, which are more likely to 
be used, should receive attention because of possible adverse effects 
on the characteristics of the longitudinal motions. 


Effect of Parameters on Stability Eoots 

klien the equations of motion are solved, the motion is obtained 
as the sum of a series of components called modes. Stability roots, 
which indicate the degree of stability of the various modes, can 
also be obtained from the equations of motion without effecting a 
complete solution of the equations. A mere detailed discussion of 
the significance of the stability roots is contained In the appendix 
of reference 4. Information obtained from the stability roots is 
most useful, when the relative magnitude or importance of the various 
modes is Known, because the roots then provide a clue to the nature 
of the complete motion. 

In the present analysis, four stability roots \ are obtained 
from the longitudinal equations of motion and are distinguished by 
the subscripts 1 to 4. The nature of the roots changes with 
variations in the parameters of the hydrofoil system. A typical 
variation in the real parts of the roots is shown in figurelO. In 
general, when the ma.gnitudes of any two real roots become equal, the 
two real roots are replaced by a conjugate pair of complex roots, 
each having the same magnitude for the real part. Thus, such pairs 
of complex roots in figure 10 are indicated by a double line and an 
3-PPropriate modification of the subscript. The magnitude of the real 
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part for such complex pairs of roots should, "b.e reaii. off the plots at 
the center of the douhle line . 

For 9vex'i’’ real root o’otninod froTii the sauaticns of motion the 
complete solution will. contain an uneriodic mode, or component, 
of tile motion. Likewise, for every pair of .complex .roots the motion 
will contain an oscillatory component. \-Then the ma^mtudc of the 
real part of any of the roots passes throvi^h i^ero, the motion hecomes 
unstahle . 

Horizontal c ente r- of - '^ravl ty location . - The effect of chanc,infT 
the horizontal location of the center of .yr'avity on the real pai'ts 
of the staMlity roots is sho-tn in figure 10 for a systera of two. , , , 
equal hydrofoils with 10 ° dihedral. For oenter-of-' gravity loc.ations 
ahead of the hydrod^.'n.amic center of the front h%drc.'foil, two reel 
roots ?]_' and and a pair of complex roots ' Xo exist. V.Hien 

the center of gravity is 2.l6c2 ahead of the front'^hydrofoii. the 

roots are unhtahle, which indicates th'^t ■'•he center of g-^avity 
has reached the forwai’d hoi’ndatr.v of the stah] e .•’egion. As the center 
of gravity is moved loa'cward. the stabili'cy slo'uly impro'ves for the 
oscillatoxy component of the motion represented ty the X--. .Jl roots. 
Meanwhile the mamaitudes of the X] and X2 roots a'Dproach each 
other and tecome equal when the center of gravity is about 1.0C]_ 
behind the front hydrofoil, vith fo.rthsi:’ resD"vard mo'^rement of the 
center of gravity the roots are coupled as tvo oscillations 
represente;i by Xt .e and then the center of gravity Is 

moved bsek to the vioinit;’^ o.f behind the front hydrofoil r.ather 

I'opid changes in coupling occur, which finally re3''.fLt in a real 
root Xj_ with a lai’ge amount of damping^ a complex pair I-2.3 uith 
moderate damping, and a real root Xl|. with slight damping; VTien tho 
center of gravity is moved back to- 'a point -5 .JiOi'jL behind the' fi'ont 
hydrofoil, the magnitude of ■ihe X;y root becomes zero, and the re.ar 
bound .ary of the stable region has been roac/ied.. 

The behavior of the loots as the horizontal location of the 
center of gravity is changed indicates that the typis cf motion caused 
by disturb<ances will be considerably inflaienced by the longitudinal 
location’ of the center of gravity. 

Pate cf change o f d o'^sTiTvm i.sh The effect on the stability roots 
of assuming the dovnwash euigle e to be 2aj^ instead of zero han be 
seen from a comparison of fi^rures 10 and 11 . No, pronounced change 
in the I'octs occurred wi-th variation in except, foi-- a shift of 

the pattern of root couplings- with respect -to, the .horizontal; center- 
of-gravity. location,* this result is consistent with indications 
obtained -earlier from a study of the in:''luence of . on- the stahle 
region. Hence, for the rest of the 'work the value .^f ■ 
assumed to be zero.' ' 


was 
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Dihedral . - The Influence on the staMlity roots of changing the 
dihedral angle from 30° to 0° is evident when figure 10 is compared 
with figure 12 . The difference in the rate at which the l 2'3 
oscillation develops with rearward center-of-gravity movement for 
the two dihedral angles accounts for the different appearance of 
the ri^it side of the diagram in the two figures. The most 
Important feature disclosed by tlie comparison is the Improvement, 
hrotight about by the use of dihedral, in damping of the component 
of motion involving the root or the complex pair 

Vertical center-of-gravity location .- Figures 13 and ik 
together with figure 10 shew the effect on the stability roots of 
varying the vertical center-of-gi-avity location from a point on a 
level with the hydrofoils to a point lOc^ above the hydrofoils. As 
had been indicated by the diagrams of the stable regions, no 
pronounced changes occur in the nature of the 2'oots when the vertical 
center-cf- gravity location is shifted. 

Rate of change of lift with immersion .- The effect on the 
stability roots of making the val\ie of §Cx/ 5 z ' twice that for 
30° dihedral is evident if figure 15 is contrasted to figxure 10 . 
Doubling the vertical-damping derivative caused marked improvement 
In the ^3.1). oscillation, which suggests that the similar improve- 
ment in damping obtained by increasing the dihedral angle from 
0 ° to 30® was a resiat of the associated increase in the value 
of ' . 


Effect of ^Parameters on Xndiclal Responses 

An indie ial ?'esPonse is the motion resulting from a unit force 
or moment suddenly applied to the hydrofoil system at zero time and 
held constant thereafter. The indlcial responses are of interest 
because they are of the same general character as the motions 
produced by types of disturbance that are likely to be encountered 
in practice . 

The longitudinal equations of motion (equations (9)) involve 
tliree variablesj hence three indlcial responses are necessary to 
define the motion caused by any specific unit disturbance. The 
three indlcial responses may be conveniently represented by the 
symbols z'^, and ©2 the change in angle of attack, 

vertical position, and angle of pitch, respectively, when the motion 
Is caused by the sudden application of a unit Cy^-force to the 
hydrofoil system. Similarly c%, z 'j^, and 0 are the response 
factors for a sudden unit Cjj^ disturbance. 
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The Indlcial responses are functions of nondimensional time s^, 
typical variations of which are shown in figtire 16. The magnitude 
of dist\irhances actually encountered, when expressed in coefficient 
form, will usually he considerably less than unity; consequently, the 
actu.al motions exuerienced will he of proportionately smaller magnitude 
than the indicia! responses presented hut will have the same type of 
variation with time. Values of the indicial responses after the 
disturbance has been absorbed by the system and new steady- state 
equ-ilibrium conditions have been reached are represented by short 
horizontal lines at the right side of the plots. Such steadj'-- state 
values are not Only new equilibrium conditions for sudden disturbances 
but also represent new trim conditions after gradual changes in the 
load condition, such as would result from the use of fuel . 

Horizontal center- of- gravlt.v locatio n.- Indicial responses for 
a unit disturbance applied to a system of two equal hydrofoils 
with zero dihedral are plotted against nondimensional time in 
figure 16 for several horizontal locations of the center of gravity. 
Values of Xj_ used in figure 16 were selected to give center-of- 
gravity locations covering all the types of root coupling shown in 
figure 12. If the center- of- gravity locations used in figure 16 
are noted on the diagram of the corresponding stable region (see 
fig. 3)^ the following points are evident: 

(1) A center-of-gravity location near the front boundary of 

the stable region is conducive to motions characterized by pronoimced 
oscillations. 

(2) A more rearward location of the center of gravity reduces 
the prominence of the oscillations but Increases the viltimate 
deviation from the attitude that existed before the disturbance. 

(3) For center-of-gravity locations near the rear boundary, no 
discernible oscillation is noted, but very large departures from the 
Initial condition occur. 

Comparison of the maximum deviations for the three center-of- 
gravity locations of figure 16 shows that, during the Interval of 
time covered by the curves, the smallest amplitude of motion of the 
hydrofoil system occurs for the case with the center of gravity back 
35 percent of the distance I between the two h3''drofolls . The 
deviation caused by a given disturbance rapidly becomes greater as 
the center of gravity is moved back of the optimum location, with the 
result that for such rearward locations a very slight disturbance 
would bring the hydrofoils to the surface or cause them to sink very 
deep into the -water. Location of the center of gravity any appreciable 
distance ahead of the optimum location appears undesirable because of 
the pronounced oscillatory motions involved. Such motions would 
be both uncomfortable and difficult to control. 
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Indicial responses for a unit Cjj^-dlstior'bance^ for the same 
conditions as for figure I6, are plotted in figure I7. The 
discussion of the effect of change in horizontal center-of- gravity 
location on the indicial responses for a unit C2*dlstur‘bance also 
applies for a unit C^j^-disturhance, with the exception that the 
amplitudes of the motions are least for the most forward center-of- 
gravity location considered^ instead of for the middle location# 

The oscillations are much more persistent^ however, for the forward 
location than for the middle location. 

Because of the large response factors involved for either type 
of disturbance, even when the best center- of- gravity location is 
selected, motions for hydrofoils with no dihedral wl.ll involve large 
amplitudes whenever a slight disturbance is encountered; hence, it 
appears evident that such a type of hydrofoil will not give satis- 
factory perf oxmance • This conclusion applies only to the arrangement 
investigated, where the hydrofoils always remain completely submerged; 
and it shoudd not be extended to cover ladder aiTangements, for which 
a change in effective area with immersion depth produces effects 
similar to those for partly Immersed hydrofoils with dihedral. 

Dihedral angle # - The effect on the indicial responses of 
Increasing the dihedral angle from 0^ to 30^ may be obtained by a 
comparison of figui’es I8 and I6 for a unit C^-disturbance, and of 
figures I9 and IJ for a unit disturbance. The figures indicate 
that the effect on the nature of the motions of changing the 
horizontal center-of- gravity location is much the same as that 
indicated in the preceding parts of the present paper. Thus, the 
most desirable center-of- gravity location appears to be about S.^Oci 
back of the front hydrofoil, as in the case for 0^ dihedral angle. 

At any particular horizontal location of the center of gravity the 
increase in dihedral causes an appreciable reduction in the indicial 
responses. The reduced sensitivity to disturbances when the dihedral 
angle was increased from QO to 30^ may have been a result of the 
corresponding increase in vertical damping. In such a case, as 
mentioned in the discussion of stable regions, a further increase in 
dihedral wo^jdd have an effect opposite to that caused by this initial 
increase in dihedral. 


pate 01 change of lift with immersion .- The effect of varying 
the rate of change of lift with immersion on the indicial responses 
for a unit Cg** disturbance may be seen from a comparison of 
figures I6, 18 , and 20. Figures I6 and 18 give the indicial responses 
for hydrofoils with dihedral angles of 0^ and 30^> respectively; 

. whereas for figure 20 the rate of change of lift with immersion 
is assumed to have a value twice that for hydrofoils of 30^ 
dihedral angle but to have other hydrofoil characteristics the 
same as for 0 ^ dihedral angle . If the case for the center of 
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gravity at 3.50o^ Is selected In each of the figures, comparison 
shows the direct relation hetween good riding character! 6tics and 
a large val\xe of SCl/Sz'. It appears, therefore, that a large 
value of ^C-j^/bz ' should he attained hy the use of arrangements 
such as hydrofoils with dihedral for which the effective area 
changes with immersion depth, or hy the use of some device that 
changes the angle of attack when the height varies. Figxare 21 
gives data corresponding to the data of figxxre 20, hut with a 
unit Cm- disturbance assumed, ResuJLts for the several center-of- 
gravity locations assumed in figures 20 and 21 Indicate the same 
influence of horizontal center-of-gravity location on the motions 
as has been shown by the computations summai’ized in figure 17* 

Longitudinal hydrofoil snacing .- Indicia! responses for either 
a unit C 2; -disturbance or a unit 0^" disturbance applied to a 
system of ’two equal hydrofoils spaced 20cq are given in figure 22. 

The horizontal center-of-gravity location in figure 22 is at 0 . 35 ^^ 
which is the same percentage of I that wa.s used in figures I8 
and 19, and other conditions are also the same as for figxires I8 
and 19. Figure 23 gives data similar to the data of figure 22 
except that the spacing has been Increased to lOOcj^. Comparison of 
figures I8, 19, 22, and 23 indicates that increasing the hydrofoil 
spacing tends to increase the restraint in pitching and thus 
reduces the response in all degrees of freedom for pitching -moment 
dist’orbances, and in all but vertical, motions for Z-foi-ce disturbances , 
The effect of increasing the hji-drofoil spacing on the motions suggests 
that the spacing should be as large as is practical, in order to 
reduce the response to a given disturbance. Figure 2 h shows the 
significance of 10c-, , 20cj^, and lOOcj^ apaclngs if the hydrofoil 
systems were att-ached to a typical flying boat. 


LATERAX MOTIONS 


Lateral stability for flying boats has not generally been a 
serious problem xxp to the present timej hence the present investiga- 
tion of the lateral characteristics of hydrofoils was brief and 
made chiefly to check the lateral stability of typical hydrofoil 
arrangements assumed in much of the study of longitudinal stabilitj'". 

In the p.resent investigation all the lateral stability 
calculations were made for a h3/'drofoil system consisting of two 
identical hydrofoils of rectangular plan form, each having 
rectanguaar tips, 30 ° dihedral, and an aspect ratio of 6. The 
center of gravity/- was assumed to have a horizontal location 0.35Z 
behind the hydrodynamic center of the front hydrofoil. The rate 
of change of downwash at the rear hydrofoil was assumed to be zero. 
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The mass of the hydrofoil system was the same as that assumed for 
the invest igation of longitudinal stahillty. The study was confined 
to what was considered the idealized case, where the supporting 
struts have no influence on the characteristics of the hydrofoil 
system. The method of treatment for the lateral motions was similar 
to that used for the longitudinal motions and is descrihed in detail 
in the appendix. 

The effects of changes in the vertical location of the center 
of gravity and changes in the longitudinal spacing of the hydrofoils 
on the lateral stability roots are indicated by the data of the 
following table : 


^1 

(chords) 

1 

(chorda) 

Lateral stability roots 

2.5 

10 

0 

-O.^hk t 1.1221 

-O.U 72 ± 0 . 20 Ui 

5.0 

10 

0 

-1.715 t 0.5231 

-0.2l;2 : 

t 0.1011 

5.0 

20 

0 

- 2 . 271 ; + 1.9581 

-0.221 

- 0.292 


The zero root that is listed for each set of values of 
and I in the table results because the system is insensitive to 
heading^ that is, the performance does not depend on the initial 
direction of travel. The remaining roots listed are either negative 
or have negative reeil parts in the case of complex roots, which 
indicates that all the systems investigated were laterally stable. 
Instability was expected in the two cases "with the higher center- 
of-gravlty location, but apparently the stabilizing effect of the 
rolling moment that is developed when the system is banked (defined 
by the value of the derivative outweighs the effect of the 

higher center-of- gravity location. Check calculations made with 
reduced to nearly zero but with other conditions the same 
as for the second case in the table showed pronounced lateral 
instability. From the foregoing results the value of 
appears to have an Important Influence on lateral stability. The 
value of this derivative is likely to depend on the depth of immersion 
of the hydrofoils; therefore it may impose a coupling between the 
longitudinal and the lateral motions and thus prevent reliable 
predictions of the lateral behavior when the longitudinal motion 
is ignored. In contrast, none of the longitudinal derivatives 
appears to be appreciably affected by lateral motions. 
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The data given in the table Indicate that raising the center 
of gravity and increasing the longitudinal spacing of the hydrofoils 
both increase the total damping in the hydrofoil system, hut the 
practical value of the increase in damping cannot he determined 
except from a study of the response factors involved. Such a study 
does not seem feasihle xintil experimental checlcs are made on the 
validity of certain of the assumptions made in developing the theory 
for lateral motions. 


SUGGESTIONS FOB FUTURE EESEABCH 


The present study is based on the assumption of small 
displacements . Because of the nonlinearity of many of the derivatives 
Involved, any appreciable departures from the assumed speed, depth of 
immersion, and other factors may cause marked, changes in the dynamic 
characteristics of the system. Studies of maneuvers, such as take- 
offs, of hydrofoil systems may consequently require step-by-step 
treatment. The development of methods of studying the combined 
motions and determination of the effects of changes in forvrard speed, 
hydrofoil loading, and moments of inertia on the motions also appears 
desirable. For seaplanes the interaction of hydrofoils, hu]-l, and 
aerodynamic surfaces must be considered. Other factors that should 
receive attention are the influence of the hydrofoil supports 
(particularly on lateral motion), the effects of power, and the 
nature of the downwash near a free surface. 


CONCLUDING BEMABKS 


A theoretical investigation was made of tandem hydrofoil 
arrangements, based on the lifting-line theory. The conclusions 
which follow apply to only the longitudinal behavior, inasmuch as 
the computations made were ins\ifflclent to .justify definite conclusions 
regarding the lateral motions . 

1 . The longitudinal hydrofoil spacing shoiild be as large as is 
feasible . 

2 . The rate of change in lift with change in depth of immersion 
of the hydrofoils should be large . Dihedral appears to be 
advantageous, if the hydrofoil is partly immersed, because with 
dihedral there is a larger rate of change of lift with change in 
immersion. The rate of change of lift with immersion will be 
insufficient for hydrofoils with no dihedral unless the area is 
composed of several panels in a multiplane arrangement. 


MCA TN No. 128‘5 


15 


3* The rear hydrofoil area should he as large as, or larger 
than, the front hydrofoil area if large variations in centor-of- 
gravity location are to he accommodated when the longitvuJinal 
hydrofoil spacing is small (of the order of 10 chords). With 
appreciably larger spacings, the arrangsTtient with the main surface 
forward appears to he sufficiently stable and should he more 
efficient than the other arrangements . 

4, The choice of horizontal center-of-gravlty location should 
he based on considerations of the resultant characteristico of the 
longitudinal motions and the hydrofoil loading. Tlie location should 
not he ahead of the hydrodynamic center of the front hydrofoil, in 
order to avoid undesirable loading. The location should he as far 
ahead of the rear boundary of the stable region as is feasible 
vrlthout incurring objectionable oeclllatlons. The best compromise 
from this latter standpoint appears to be a location near the center 
of the stable region. For tvro equal hydrofoils in tandem the best 
location appears to be back aboiit 35 percent of the distance between 
the hydrofoils. 

5* If the effects of power are neglected, the vertical centcr- 
of-gravity location appears to be of little importance, low locations 
being somewhat advantageous, 

6. A reduction in tho pitching radius of gyration will cause 
an approclab3.e increase in the range of horizontal center-of-gravlty 
location that will bo stable. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., May 9, 19^7 
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APPSWDIX 

METHODS OF THEORETICAL TREATMENT 
SYMBOLS 

h.c. 

hydrodynamic center 

c.g. 

center of gravity of hydrofoil system 

X— , Y-, Zr-axes 

rectilinear reference axes fixed in hydrofoil 
system, with origin located at center of 
gravity (The X-axis is alined in the direction 
of the initially undisturbed motion. The Initial 
position of the Y-axis is directed horizontally 
to the right. The Z-axia is directed do\mward. ) 

X, Y, Z 

forces along X—, Y-, and Z-axes, respectively 

L, M, W 

moments about X-, Y-, and Z-axes, respectively 

Z'-axis 

axis, directed vertically downward with respect 
to the earth from origin located at center of 
gravity of hydrofoil system 

Z* 

displacement along Z'— axis 

CD 

angular dlsplacments of reference axes from 
initial positions, radians (see fig. 25) 

a, 3 

angles, in radians, giving instantaneous orientation 
of reference axes with respect to path of motion 
(see fig. 25)ji thus a is angle of attack 
and 3 angle of sideslip at center of gravity 

V 

linear velocity of center of gravity 

n 

angular velocity of hydrofoil system about center 
of gravity, radians per second 

U, T, w 

components of V along X—, Y-, and Z-axes, 
respectively 

Vi <1> 

components of Q about X—, Y-, and Z-axes, 
respectively 
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W w©j.ght of hydrofoil 33’-stem 

m mass of h^'^drofoll system 


Pw 

1 

2 

S 






6 


^<1 

3n 


radii of gyration of hydrofoil system about 
respective reference axes 

density of v^ater 

subscript used to designate front hydrofoil In a 
system of two hj'-drofoils in tandem 

stibscrlpt used to designate rear hydrofoil In a 
system of two hydrofoils in tandem 

total projected area of immersed part of hydrofoil 
system under conditions of stead;' undisturbed 
motion 

total projected area of nth hydrofoil 

chord of nth hydrofoil 

span of nth hydrofoil 

aspect ratio of nth hydrofoil 

dihedral angle of nth hydrofoil, in radians unless 
specified othen-rlse 

dihedral- angle when angle is same for all hydrofoils 
in system 

angle of attack at hydrodj'namic center of nth 
hydrofoil, radians 

Induced angle of attack at hydrod' mawi-f c center 
of front hydrofoil, radians 

downwash angle at hydrodynamic center of rear 
hydrofoil, radians 

rate of change of e with a. 

rate of change of $ with qc-^/V 

angle of sides.lip at hydrodynamic center of nth 
hydrofoil, radians 
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nondimenslonal rolling Telocity at hydrodynamic 
center of nth hydrofoil, based on local rolling 
velocity in radians per second, hj^, and T 


nondimensional yawing velocity, with definition 
similar to that for 


( 1 ) 


lift on hydrofoil system, measured at center of 
gravity in direction pei*pendicular to V and 
converted to coefficient form hy dividing 


lift on nth hydjTofoil, measured at hydrodynamic 
center of hydrofoil under consideration in 
a direction pai’allel to C-r and converted to 

1j 


coefficient form hy dividing hy 


2 \v 


lift on nth hydrofoil, measured at hj’’dro dynamic 
center of hydrofoil under consideration in 
direction perpendicular to local relative motion 
and converted to coefficient form hy dividing 



drag on nth hyrdrofoil, measxn'ed at hydrodynamic 
center of hydrofoil under consideration in 
direction parallel to local relative motion 
and converted to coefficient form hy dividing 

weight of hydrofoil system converted, to coefficient 


form hy dividing hy 


Ip V^S 
2 w 


side force on hydrofoil system, measured at 
center of gravity in direction of Y-axis and 
converted to coefficient form by dividing 


^7 


side force on nth hydrofoil, measured at hydro— 
dyns.^tuic center of hydrofoil under consideration 
in direction parallel to Y-axis and converted 

to coefficient form hy dividing hy g^w^^^n 


coefficient of Zr-force^ vith definition similar 
to that for Cy 

rolling moment about Z-axis, converted to coefficient 


rolling moment at hydrodynamic center of nth hydrofoil 
about axis parallel to X— axis, converted to 


pitching moment aboiit Y-axis, converted to coefficient 
form by dividing by 

coefficients of yax^ring moment, with definitions 

similar to tliose for Cj and , respectively 

' 'Ti 


X— component of distance from center of gravity to 
hydrodynamic center of front hydrofoil, Ct_— units 

X— component of distance from hydrodjTiamic center of 
rear hydrofoil to center of gravity, Cy-unlts 

distance between hydrodynamic centers of the two 
hydrofoils measured paa’allel to X-axis, Cj|_— units 

Z— component of distance from center of gravity to 
hyarody::ariic center of nth hj^drofoil, Cj^— luiits 

operating depth; distance from water surface to 
hydrodj^i.-ntiic center of nth hydrofoil during 
steady undisturbed motion, Cj^— units 

operating depth when depth is same for all hydrofoils 
in system 

parameter of nth hydrofoil used to determine value 


form by dividing by ip^V^Sb^ 


coefficient form by dividing by 


the derivative 



of 
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Y-conponent of dlstsncc from hydrodynamic center 

to centroid of lift on one panel of nth hydrofoil, 
'bn"'^ts 

vertical displacement of center of gravity during 
disturbed motions, C]^ -units 

vertical displacement of hydrodynamic center of nth 
hydrofoil during disturbed motions, Cj^-units 

imss of hydrofoil system, -|p^SC]_ -units 

mass of hydrofoil system, yp^b^-units 

radius of gyration about Y-axis, c^^-mlts 

radii of gyration about X- and Z-axes, respectively, 
b^-vuaits 

time, seconds 

time, c^/v units (To convert nondimensional 
time into second units use t = SgC^^/V. The 
time scale may alternatively bo converted into 
distance traversed if values of Sp are multiplied 
ty ciO 

time, bj^/v units 

stability root, with various mmnerical subscripts 
used to distinguish the different roots 

disturbance function; a Z-force of variable magnitude, 
time history of which is Indicated by form of 
f\mction (The complete description of any arbitrary 
disturbance acting on the hydrofoil system may be 
expressed by use of this and the additional 
disturbance functions M(t), Y(t), L(b), and 

N(t), with definitions similar to that for Z(t).) 

nondimensional disturbance function, similar 

to Z(t) but with force expressed in coefficient 
form and with time in nondimensional units 
(Similar definitions apply to Cjh(Sq), Cy^s-b), 

Cz(s-b), and Cjst,).) 
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indlclal responsec giving motions in a, z* , 
and 6, respectively, caused ty sudden 
application of unit C^-disturtance to 
hydrofoil system 

Indlclal responses giving motions in a, z’, .and 9, 
respective.ly, caused by sudden application of 
unit Cjjj-dlsttu'bance to hydrofoil system 

empirical constant used to determine value 

empirical constants used to determine value 
of 

empirical constants used to determine value 
Of S(CB)^/aa„ 


Longitudinal Equations of Motion 

The longitudinal motions of the hydrofoil system are referred 
to the system of axes described in the list of symbols. The choice 
of axes that correspond to those customarily employed in studies of 
airplane stability should facilitate extension of the present hydro- 
foil theory to include the effects of aerodynamic surfaces. The 
equations of motion are based on the assumption that the hydrofoil 
system can be replaced by a particle at its center of gravity having 

, about the respective 

reference axes equal to those of the hydrofoil system. The analysis 
is also based on the assumption that the velocities V in the 
direction of motion and u along the X-axis are constant and that 
departures from the initial conditions of motion eire small. The 
further assumption is made that the longitudinal displacements Z’, 

9, and along the Zr-axis, in the plane of symmetry of the hydrofoil 
system, are Independent of the lateral motions involving the 
displacements 0, t, and along the Y-axis. This assumption yields 
satisfactory theoretical predictions of the motions of airplanes in 
normal flight and appears warranted, based on the nature of the 
deviations involved, in the treatment of the longitudinal motion 
of hydrofoils. Its application to the lateral motions of hydrofoils 
is made with reservations, as mentioned in the main test. 


a mass m and radii of gyration k^^-, k^ 




<%, 0 . 


m 


k. 


^ 2 ^ ^3 


ki|., k^ 
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By the use of D 'Alemherts principle, the f ollo^rinc equations 
of equilihriuQ at the center of gravity are inritten for the forces 
and moments involved in the longitudinal motions: 


d‘ 

a— 

dt' 


- - mqV = + Z -I- + qS + Z(t) 

2 ^ Sw !^iZ ' ee dq 




ink- 


, 45 * z 'ciL + 9 ^ H M(t) 


dZ ’ c9 


( 1 ) 


where Z(t) -and M(t) are arhitrai'y distvirhance functions. The 
equations have the same form as the familiar equations of longitudinal 
motion for an airplane, except for the addition of d.erivatives with 
respect to Z ' and 9 . The equation of equi.lihi-iimi involving the 
X-forces is oiaitted hecause is assumed constant. Equations (1) 

can he simplified hy using w = M, a = -, and q = — to give 

dt V dt 


..aa AQ 


CGcr— -5- 
ca 



^ (jZ d9 5Z 

'’3? dt s; 




-I- 

da 


Z'M_ 

^Z ' 


d0 

d9 ^ dt dq 


M(t) 


( 2 ) 


If equations (2) are reiiTitten in a nondimens ional form, the soluticns 
obtained vdll he general in character. The method used to make the 
various terms of the equations nondimsnsional involves expressing all 
angles in radians, all forces and moments in the standard HACA 
coefficient forms 
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M 




( 4 ) 


all lengths in terms of the chord of the front hydrofoil, all 
times in terms of the time c-j^/V required for the system to traverse 

the distance c^ along the path of motion, and the mass in terms 

of •units. The nondimenslonal quantities of mass 

time Sq, vertical displacement z*, and radius of gyration Ky 
about the Y-axis thus hear the folio-wing relations to the 
corresponding dimensional quantities ; 




-EL 




" oi/y 




Ky = 

^ Cl 


(5) 


(6) 


(V) 


( 8 ) 


In equations (3) to (8), is the density of -water and S is 

the total projected hydrofoil area in the hydrofoil system. 

The nondimenslonal form of equations (2) becomes • 




da 


d0 


^sc d.S(j 


da 


+ 0^ + ^ _Z + c (s.) 

d0 ds^ Z' c^ 


V 




?d^9 

^c —“2 
ds 


Ufi Ey 


dC. 


m 


+ z*^ + 0 ^ + ^ 5PeL_ 


da 


da* 


d9 


dsc .q-Cn 


Cm(sc) 


( 9 ) 


2h 
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Also, from geometric considere.tions, 


dBc 



( 10 ) 


In equations (9)., functions of nondinensional 

time that describe the application of disturbing force and moment 
coefficients to the hydrofoil system. The methods used to make the 
terms of equations (9) nondimenGioneJ. have the advantage that the 
nondimonsional equations obtained retain the same form as the original 
force equations) consequently the physical significance of the non- 
dimensional equations should be more readily evident. Solutions of 
motion obtained from equations (9) are llke'/ise nondimensional and 
may be considered as proportions, applxcablo to all similar hydrofoil 
systems, and capable of conversion to customary engineering imlts in 
any given case bjr use of the characteristic dimensions c^ and V 
pertinent to the specific design. 

Stable regions and stability roots for the longitudinal motions 
can be obtained from equations (9) in con^Juncticn -with equation (10) 
by methods discussed in reference k. The stability equation for the 
longitudinal motions has the form 


Boundaries for the stable regions vere obtained from the conditions 


for the divergence boundary. The quantities involved in equatloias (12) 
and (13) are the coefficients of equation (ll), which in tur-n are 
functions of the factors of equations (9) and (10). G?hus, 


aD^ + bD^ + cD^ + dD + 0 = 0 


( 11 ) 


(be - ad)d - b^e = 0 


( 12 ) 


for the oscillatory boundary and 


e = 0 


(13) 



( 14 ) 
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Equation ( 12 ) Is the familiar Rauth*e discriminant, but Its expression 
In terms of the factors In equations (9) and (10) is considered too 
lengthy to be presented here. 


Longitudinal Derivatives 

Values must be assigned to the varioxis partial derivatives 
appearing in equations (9) before the equations can be solved. No 
experimental values for the derivatives were available^ hence computed 
values were used. The computed derivatives were eTOluated on the 
basis of experimental hydrofoil data obtained from restate of tests 
made in the Langley tarJc no. 1 at various immersions and speeds. 

A discussion of the methods used to compute the various derivatives 
follows. Data presented in connection with the discussion are for 
hydrofoils of rectangiaar plan form and tips, with an aspect 
ratio of 6, and operating at a velocity of 20 feet per second. 
Experimental restate indicate tliat, for a given angle of attack, 
marked changes in the lift and drag coefficients of hydrofoils 
occur with changes In speed. The values of the derivatives would 
undoubtedly be equally affected by any pronounced change in speed 
from that assumed in the investigation. 

Chang e in ^fo rc e with v er tlcal displacement of the center of 
sravity .~ If the center of gravity moves downward, the 

hydrofoils are immersed deeper in the water. Experimental restats 
indicate that an increase in the depth of Immersion of a hydrofoil 
is accompanied by an increase in the magnitude of the lift obtained. 
The increase In lift is proportional to, and of the same sign as, 
the initial lift. Thus, 





(15) 


Values of k^ are given in figure 26 for a dihedral angle of 0 ^ 

and in figure 2 ? for dihedral angles of 20° and 30°. The value 
of kj^ depends on the normal operating depth z^^ of the hydrofoil. 

The discontinuities in the curves of figure 2 ? coincide with the 
point where the tips of the hydrofoil break the surface. Hx fig- 
ure 26 and subsequent figures (Cl) is based on the total area 

of the hydrofoil instead of the immersed area, and z« is 

meaexzred in chord lengths of the particular hydrofoil under 
consideration. 
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The Talue of fo- a coiaplete hydrofoil system is the 

negative sum of the values of oC^^^dz* for the individual hydrofoils. 

The values of 8Cj^y/8z» for the various hydrofoils are derived from 

the ■’■allies obtained from figures 26 or 2? hy making 

proper allowance for the different areas and chords that are used 
to males the vaj’ioui terms nondimensional. 


Change In Z-force with angle of attack The value of 

the derivative SC^/^ct is the negative sum of txie values 
of (that is, the elopes of the lift curves) for the 

individual hydrofoils. As in the case of ^Cg/8z’, differences 
in the areas used in forming the coefficients must he taken into 
account when the addition is m£ide. The slope of the lift curve 
depends on the depth of immersion of the hydrofoil. Typical 
variations of the’ slope are given in figure 28 for 0° dihedral 
angle and in figui'e 29 for various dihedral angles. When figures 28 
and 29 a.re used to determine the slope of the lift curve for the 
rear hydrofoil, the value obtained is with respect to the local 
angle of attack etp at the rear hydi’cfoil. In general the value 
of Op is less than that of a (measured at the center of gravity) 

hy the amount of the downwash angle € at the rear hydrofoil. The 
slope of the lift curve for the rear hydrofoil must he corrected 
for down^;ash to give the required slope with respect to a. The 
correction is applied hy multiplying the slope obtained from 
figure 28 or 29 hy the factor 1 - where has some value 

in the range 





Ai 


(16) 


In equation (I6), A3_ is the aspect ratio of the front hydrofoil 

end 8(C^) l^o^i is the lift-curve slope obtained from figure 28 

or 29 for the front hydrofoil. 

Change in Z-force with pitch attitude A change in 

the pitch"^HTude of the hydrofoil system will cause a differential 
change in the depth of immersion of the hydrofoils. The effect on 
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the Z-force may he estimated from the geometry of the system and 
the data of figures 26 and 27) thvis, for two tandem hydrofoils 




Sp 

Xt ^ 

1 .q 


1 

2 8z'o 


X2 


(IT) 


Chnnge in Z-force with pitching ve]-Ocfty The chief 

effect " of ""aTpitchlng velocity ahout the center of gravity of the 
hydrofoil system is to cause a change in local angle of attack at 
each hydrofoil. The change in effective camber for the pitching 
hydrofoil introduces a small additional component of vertical force, 
(See reference 5 -) The total effect for two hydrofoils in tandem 
may be assumed to be 




( 18 ) 


Where 




(— + 0.5) 


(19) 



In equations (18) and ( 19 )^ j ^2 lif't>-curve slope for 

the rear hydrofoil, based on the local angle of attack 02j and 


% 


28 NACA TN No. 1285 


Xg are the X-conrponents of the locations of the front and rear 

hydrofoil hydrodynamic centers from the center of gravity expressed 
in terms of Cji and e indicates the rate of change of downvaeh 

angle at the rear hydrofoil with change in nondlmensional pitching 
velocity qC]_/V. The value of 6 ^^ will he in the range 


(a, 

^ jtA^ 

Chang e in pitching mome nt w ith ve rt ical d l sni .ace ment of the 
center~of gi’avity The changes in lift, mentioned in the 

discussion of the change in Z~force with vertical displacement of 
the center of gravity, produce moment changes about the center of 
gravity, the magnitude of \ 7 hich depend on the X-components of the 
distances of the hydrofoil hydrodynamic centers from the center of 
gravity. The drag also increases with deeper immersion of the 
hydrofoils. Analysis of the data obtained in Langley tanlc no. 1 
indicates that the change in di’ag can be expressed as 



n 


Values of ^ 


3 


are given in figure 3O 0 ° dihedral 


angle and in figure 3I for 30° dihedral angle. The drag changes 
multiplied by the Z-ccmpcnents of the distances from the center 
of gravity to the hydrofoil hydrodynamic centers give the drag 
contributions to the change in pitching moment. For two hydro- 
foils in tandem 


5c 

m 

ozf 


S 


^2 Cl 




'2 


Xg- 


{%) 




^ fa)2 


=>2 c . 


S ^2 8 z» 


(23) 
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Change in p i tching moment with angle of attac k bC^/ba. — 

Physical considerations lead to the expression, for two hydrofoils, 


^5 

8a 


where 


*=> oa-j^ 


, ^(^1)2 

- ea) . 2:2 + 

dOg 

So , 

+ _2 (1 _ .) 

CVI 

^(^d)2 

ba^ 

(Cp)j^ is the drag coefficient of the 
area' of the front hydrofoil.; ^Cp'^ 


S 


rn ) _ 


( 24 ) 


Of the rear hydrofoil hased on the area of the rear hydrofoilj 

and Zj and Zq are the Z-components of the locations of the front- 

and rear-hydrofoil hydrod’/namic centers from the center of gravity, 
expressed in terns of cj. 

The slope of the drag curve for each hydrofoil must he Icno’Tn to 
determine bCjj^/ba from eqi;ation ( 2 k). The empirical relation 



(25) 


was obtained from an analysis of the experimental data. Values 
of kij^ and k^ varied with the depth of immersion of the 

hydrofoils in the manner shown in figur-e 32 for 0° dihedral angle 
and in figure 33 for 30° dihedral angle. 

Change in pitching moment with pitc h attitude bC^bO.- The 

differential change in the depth of immersion of the hydrofoils 
introduced hy a change in the pitch attitude of the hydrofoil 
system leads to variations in the lift and drag for each hydrofoil. 
These variations can he translated into a variation in pitching 
moment about the center of gravity hy use of the geometry of the 
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hydrofoil system and equations (l5) and (22). For two tandem 
hydrofoils 


52m 

hd 



!g !i ^ ^ X,., 

S Cg ^ ^ ^ ^ 


Sg 

s 


’^ 2^2 


(26) 


I 

Change in pitching moment with pitching velocity /5-4.- The 

oELly important contritution to the pitching ’moment^bduced by a 
pitching ve?-Ocity about the center of gravity is that associated 
with the change in lift on each h^-arofoil as a result of the change 
in local angle of attach. Thus, 


_ £1 X a 


V 




(27) 


Lateral Equations of Motion 

Equations expressing the equilibriums of the forces and momenta 
involved in the lateral motions are trritten on the same assumptions 
as those used to obtain the ].ongitudlnal equations. The equations 
of lateral motion are 
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+ mrV = ^ (— + 

dt2 St J 


+ r|^ + Y(t) 
8p or 


rak^2 d^ = + p^:; ^ ^ 

dt^ ^ S0 Sp Sr 


(28) 


mk^ — * = V — + + p^ + r — + H(t) 

^ dt^ 8v ^ Sp 


vhere Y(t), L(t), and N(t) are arbitrary'- disturbance ftmctione. 

Equations (28) can be simplified bj using v = — . 0=^, n=.^ 

^ ^ dt V ^ dt 


and r = ~ to give 
dt 




mvffi + mv52 = p.5^ .K,) /k + w ) + M + a 2K + Yl-t) 

dt dt ”313 U 7 * at ap dt Sr ^ 




dt' 


S0 


dt 8 p dt dr 


= „/t) 

^ dt^ S0 ajiJ at dp dt dr ^ ^ 


(29) 


Equations (29) will next be written in a nondlmensional form similar 
to that used for the longitudinal equations. Thus, al l angles will 
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te expressed in radians and all forces and moments in the standard 
NACA coefficient forms 


Oy = — S— 

i<V'f2s 






c, = L. 




(30) 


(31) 


(32) 


n = 1L__ 

^T=0bi 


(33) 


Because of the different basis for forming the moment coefficients 
(cf. ecjuation (4)) in the nondimensional lateral equations of 
motion, all lengths will be expressed in terms of the span of the 
front hydrofoil bj^, all values of time in terms of the time b^/V 

required for the system to traverse the distance b^ along the 
path of motion, and the mass in terms of units. The 

nondimensional mass u-^^, time s-]^, and radii of gyration Kjj- 
and thus boar the following relations to the corresponding 
dimensional quantities: 


M-t = 




( 34 ) 
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t 

^l/Y 


(35) 



(36) 



(37) 


The nondimeneional form of equations (29) "becomes 


W, «%/ 




V 


MTd%' 


2 + M. 5i!_ + 2L + c,{8t) 

j__ 2 ^ ^rt( ,i« Tihn iio, rb", i b 






Y 


f^b^z' 


2 £i_ 

dspS 


0 ^ + o& + ^ p (^) 


where Cy(sb), Cj(a^), and Cj^(sp) are functions of nondimensional 

time that can be used to define the application of any lateral 
disturbance to the hydrofoil system. 
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Lateral Derivatives 


In order to obtain a solution from equations ( 38 ), the various 
partial derivatives involved must he given numerical values. Wo 
experimentally determined values were available for an;’’ of the 
derivatives^ and computed values were therefore used. Experience 
has shown that theoretical methods are lun’eliable for obtaining many 
of the lateral stability derivatives of airplanes. This fact, 
coupled with the additional complication of the presence of a free 
surface, suggests that theoretical computations of the derivatives 
for hydrofoils will be even less satisfactory. Elaborate theoretical 
analyses to obtain the "values of the lateral stability derivatives 
of hydrofoils, therefore, appear to be unjustified until experimental 
data are available for use in checking the accuracy of computed values. 

For most of the lateral derivatives, the values of the 
derivatives were first computed with respect to the hydrodynamic 
center of the hydrofoil for motions at the hydrodynamic center j from 
the geometry of the hydrofoil system ’biie der-?vatives at the center of 
gravity of the hydrofoil system for motions at the center of gravity 
were obtained. The following discussion will be mainly confined to 
methods of computing the lateral derivatives at the hydrodynamic center 
of the hydrofoil. Such derivatives can be readily converted to 
derivatives at the center of gravity of the hydrofoil system by the 
use of elementary mechanics when the geome’try of the system is known. 
Wmerical data presented in connection with the discussion of the 
lateral derivatives were obtained from the same sources and the same 
operating conditions as those used in obtaining the longitudinal 
derivatives. The expressions derived are for the lateral derivatives 
of an "ideal" hydrofoil system without sLxpporting struts. The presence 
of the supporting struts usually required will imdoubtedly have a 
large influence on the values of certain of the lateral derivatives. 

Change in Y— force with sideslip c)0y/S3.— During sideslip the 

effective angle of attack is differentially altered on each side of 
the hydrofoil, which changes the lift on each half in such a way that 
a component of side force is introduced. This effect is a function 
of the dihedral of the hydrofoil. In addition, the direction of the 
drag force is rotated to one side during sideslipping. The sum of 
these effects is 



( 39 ) 
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where 


coefficient, based on of the Y-coraponent 


n 


of force at the hydrodynamic center of the nth hydrofoil and 3. 

is the sideslip angle at the samo point. The dihedral angle of 
the nth hydrofoil in radians is indicated by The value of 

^ (pL)n 

or 29, and the value of (Pd)j^ given in figure 3^ for 0° dihedral 
angle and in figure 35 for 30° dihedral angle. 


required in equation ( 39 ) can be obtained from figures 28 


Change in Y-force with angle of b ank 8 Cy/ ^ 0»~ The value of the 
derivative was estimated by treating each panel separately 

as a hydrofoil of which the dihedral angle, angle of attack, centroid 
of lift, lift— curve slope, and immersed area vary with angle of bank. 
The change in effective aspect ratio, which should be small fer emal 1 
changes in banlr angle, was neglected. The variation In dihedral angle 
and Immersed area with angle of bank was obtained, by graphical 
methods, for banking about the center of gravity of the hydrofoil 
system. The changes in lift-curve slope and centroid of lift with / 
dihedral angle were obtained from figure 36. The value of 730^ 

in this figure is for a lift coefficient based on the projected area 
of the hydrofoil while banked, rather than on the initial projected 
area, and vrith the lift measured vertically regardless of the bank 
attitude. The lateral displacement of the controid of lift from the 
Juncture of the hydrofoil panels is given by the value of y^,^ in 

figure 36. In order to make y^,^ nondimenslonal it is expressed in 

terms of twice the projected span of the banked panel. The new angle 
of attack of the panel after a change in bank is 


a = Oq cos Fq sec F (40) 


where the subscript 0 refers to the initial values for the hydrofoil 
panel, and f and a are the values of the dihedral angle and angle 
of attack of the panel after a change in bank. (Note that F = Fq ± 
where the sign depends on whether the left or right panel is involved. ) 
Eq.mtion (40) and the values of and y^,^ obtained from 

figure 36 can be used to determine the magnitude and point of 
application of for each banked panel. The value of Cy for 

the banked hydrofoil is then determined by rules of simple mechanics. 
Tlie value of SCy/S^ is obtained graphically by plotting the values 

of Cy determined for several values of 0 and measuring the slope 
of the resulting curve. 
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Change in Y-force with rolling Telocity ^ estimation 

.... — ,, — . — ■ — -«* — 

of the value of the derivative was obtained on the 


assumption that the side force would he zero for rolling of the 
hydrofoil about its effective center of curvature in front elevation. 
The derivative for rolling about the center section of the hydrofoil 
can then be obtained by an expression of the form 


<t) ‘ ' 




(hi) 


The parameter rj^ is given in figure 37 for various dihedral angles. 


rbi 


Change in Y— force with yawing velocity SCy/S-:^.— The derivative 
rb-i 

^CY/h-— was assumed to be zero for yawing about the hydrodynamic 
center of the hydrofoil. 


Change in rolling moment with sideslip SC^/SP.— The differential 

change in lift, produced on each panel of a hydrofoil during sideslip, 
introduces a component of rolling moment about the center section. 

An additional component of rolling moment arises because. the point of 
application of the side force produced by sideslip lies above the 
center section. The sum of these effects is 




n 


^ n ^ ^ 



(te) 


where yf. 

'^n 


is obtained from figure 36 and 


from figure 29. 


Change in rolling moment with angle of bank 8 C^/^ 0 .— Increments 
of Cj^ and Cy, ""caused by a change in’angle of" banSj "can be computed 

by methods outlined in the discussion of SCy/^ 0 . These increments, 
when multiplied by appropriate moment arms (expressed in span 
lengths), are used to obtain a plot of C^ against 0, from which 
the value of SC^/S^ is measured. 
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Change in rolling moment with rolling velocity 
Eeference 6 gives -0.2 as an average value of the derivative 



for a conventional airplane wing. 


The value for a hydro- 


foil will probably he somewhat smaller, hut in the absence of 
experimental data the average value mentioned was used for rolling 
of the hydrofoil about its center section. 


Change in rolling moment with yawing velocity 
averago"value " 


aoj/aSh.- The 


■ fa)n 

n 


(U3) 


was used for the derivative 


30,/arb. 


Reference 6 Indicates that 


this value is suitable for wings with moderate taper, and the loss 
of lift on parts of a hydrofoil that approach the surface would 
result in a similar lift dlstrihution if the hydrofoil had dihedral. 


Change in yawing moment with sidesli p SC^/Sp.- During sideslip 

the lift vector for each panel of a hydrofoil remains perpendicular 
to both the hydrofoil leading edge and the direction of motion. 
Hence, the projection of the lift vector on the horizontal plane " 
rotates forward for the leading panel and rearv/ard for the trailing 
panel. The resulting couple about the hydrodynamic center of the 
hydrofoil is 



n 



tan 




Change in yawing moment with angle of bank If, during 

hanked motion of a hydrofoil, the centroid of drag fOT~each panel is 
assumed to have the same location as the centroid of lift and if the 
additional assumption is made that the variation of drag with lift is 
the same in the banked attitude as for zero banlc, can be 

confuted by methods similar to those used for BCy/S^ and oC^/^^. 
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Change in yawing moment with rolling velocity The 




average value given in reference 6 for an elliptical distribution 

of lift was used for the derivative SC Thus 

7 Y 



(45) 


The elliptical loading was assumed to approximate the loss In lift 
over the tip parts of a hydrofoil with dihedral and with the tips 
at the water stirface. 


Change in yawing moment with ya’-rlng velocity 

value 



The 



(46) 


appears to he a suitable approximation to the expression given by 
Glauert for elliptical v?ings (see reference 6) and hence was used in 
the calculations. The selection of elliptical loading was based on 


the same considerations as for the derivative 
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Z^/yc/re Z.- Definiiion of symbo/s for a represenfa- 
t/ve hydrofoil system. 
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Figui^e 3- Stahls regions for three area disirtbufion^. r =0°;e*0; 
X-IO.Oc,\ z^=/. OOc^; Ky = 6.67Cf. (For th/s f/gare c, /•$ the 

chord for arrangement /j //? /vh/ch the f/vo hya/rofo//s are 
/denf/ca/.) 


NACA TN No. 1285 Fig. 


N 



n 


Fi^'jre,^- ^tabJe re^/ons for two doconcoast) an^i^s - 
I = 1 0.0c/, Zq = /.OO c, • t<Y = 6 S7 c, * 


= >S 


/ J 


NACA TN No. 1285 



Future 5~ 


Sto^/e 

r 


MUVI^UKT 

COMMITTEE FO# AERONAUTICS 


reg/ons For three, /m mers/or? c/epths 
e = ^; I = /O.O c, ^ Ky = 6.67c, 




cn 


NACA TN No. 1285 Fig. 



y NATIONAL ADVISORY 

^ ' COMMITTEE FOft AERONAUTICS 
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Figure 7- SfodJe r'&g/on wit/? 1Ci/'h2.' doub/e fbot for 
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Figure 10.- VartQf/on of rcc/t parts of 
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Figuf^e 12.- Vor/at/on of reo/ parfs of 

5 fat>i//ty roofo cu/fh X, . n=0° ^ 

•5^ — ) € — O f / — /O-O Cy 'j — hOO C/ j 

Ky = 6.67 c, 3 - 


Magn/tt/d£ of^ reo/ part^ of 3tat>///fy roo-f^ 
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3 to b/ /tty roots co/t/> y, . 30°^ 
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3 tob///fy roots cu/th X/ • - 30° j 
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Figure 15.- Far/of/on of re a/ par of 

sfab/Z/fy roofs cu/ff) x, for /n c r eased 
vert/caf damping . S 2 - S, \ € - O 

I - ZO.O c, /<Y = 6^G7c, , = S. OO . 
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F/gure /&.- Inc/zc/o/ respo/7ses fo, un/t d/sfo^bonci 2 . n=Oj 

e^Oj l^/O.Oc, ; Zg^/.OOc,) Ky^6.67c,-z,--5.00cp^ x, = ~/.Z5c^_, 
3.50c^j or4-.80c^. 
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F/gure n— Inc^/c/oZ responses for un/Z dtsf urboncQ . F- 0 

S^ = S,;6=0j 2=/0.0c,j z^=WOc,i /<y=6.67c,;, z,= 5.00c,- X, = ~l^5c,j 

3.50c, j or 4.80c, . 
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f-igu/'B !3.~ f/7(//c/o/ /"espo/yses fof un/t p/3turi>ance . n-SOj = j 
e ^Oj l^lO.Oc,^ T^=n^cpj Ky^6.67c,j zp3.00c^j x, =-2^.00 
S.SOc,^ 4.50 c, j or 5.25 c, . 
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Figure /9.-~ Jndicio/ responses for umf c/isfurbance. 

^=Oj l^/O.OCf) 7L^^l.l4Cfl Ky^&61Cfj -L^^^.OOOfj -2.00c^y 
3.F0c/ y 4.SOc^ y or F.Z5c^> 
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Figur'e2.0.~ Incfic/o/ responses for un/f disfurbonce. 

doub/e +bat ior P^30°j 5 ^= 5 ,; e^O; ? = /O.Oc,J Ay=d:d7c,; 
■z^Js.OOc,) X, = -Z.OOc^j 0.80 c,^ 3.50c,, or 6.Z5c,. 
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Figure 2/L— Incf/r/oZ respo/?ses for un if C/ryj d/stur bonce . 

ICiy'hz double fbot for r^JO° e^Oj l^lO.Oc,j Ky ^6.67Ci‘j 

< - -Z.OOc,, 0.60c,j 3,50Cf^ or 


z^ = d.OOc,^ 


^or 2- 
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Figure ^Z.- Ind/c/aZ responses for unit C 2 _ond umf d/sfurbonces. 

r=3dj e = 0- l^ZO.Oc,: 7-0= n^c,: ky=&67c,-^ -Z.,^5.00cp^ 

= 1.00 Of . 
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Figufe Z3 - Ind/c/o/ responses for ur>it C 2 ^ and umf C/-r> disfurbor>ces. 

r=30° S^ = S,j €=Oj l=/OO.Oc,: 1.74c, : Ky = 6.67c,; 00c;, 

= 33,00 . 
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F/^t/re Posit n/e iS(S/7sos of axe>s and motions. 
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